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ABSTRACT 

Nonlinear algebraic functional expansions are used 
to create a form for the unsteady aerodynamic response 
that is consistent with solutions of the time-dependent 
Navier-Stokes equations. An enumeration of means of 
invalidating Frechet differentiability of the aerody- 
namic response, one of which is aerodynamic bifurcation, 
is proposed as a way of classifying steady and unsteady 
aerodynamic phenomena that are important in flight 
dynamics applications. Accommodating bifurcation phe- 
nomena involving time-dependent equilibrium states 
within a mathematical model of the aerodynamic response 
raises an issue of memory effects that becomes more 
important with each successive bifurcation. 

NOMENCLATURE 

C L lift coefficient 

C t rolling -moment coefficient 

% Laplace transform 

p surface pressure 

s Laplace transform variable 

Re Reynolds number 

t observation time 

U 0 axial velocity at wing center of gravity 

u,v fluid velocity components In x,y directions 


^Staff Scientist. 
Associate Professor. 


v Q vertical velocity at wing center of gravity 
X,Y space-fixed coordinates (Fig. 1) 
x,y moving coordinates (Fig. 1) 

x,r,a moving coordinates (Fig. 3) 

a angle of attack (Fig. 1) 

v kinematic viscosity 

E running variable in time 

x time at origin of step 

* stream function 

roll angle 
u vorticity 

INTRODUCTION 

There is now widespread recognition that the evolu- 
tion of modern aircraft design makes mandatory an ever 
closer interaction between workers in the fields of 
flight dynamics and unsteady aerodynamics. Testimony to 
this recognition was the recent Joint sponsorship by the 
AGARD Fluid Dynamics and Flight Dynamics Panels of a 
meeting entitled "Symposium on Unsteady Aerodynamics- 
Fundamentals and Applications to Aircraft Dynamics" 

(J[) . (See also the useful comments and recommendations 
by the authors of an evaluation report of this symposium 
(2)). In a contribution - to the symposium, Tobak and 
Chapman (J) continued their effort to provide a consis- 
tent formulation and theoretical method for studying 
nonlinear problems in flight dynamics. They focused 
attention on the role of bifurcation theory in the math- 
ematical modeling of the aerodynamic contribution to the 
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aircraft's equations of notion. Demonstrating the 
necessity of accomodatlifg aerodynamic bifurcation phe- 
nomena within the mathematical model raised a number of 
themes that touch on the concerns of workers in both 
flight dynamics and unsteady aerodynamics. In the 
interest of widening the ground for common endeavor, an 
elaboration of three of these themes will be the subject 
of this report. 

First, we further explore the idea of using nonlin- 
ear algebraic functional expansions as a more concrete 
means of representation than the abstract use of func- 
tionals themselves. We have already noted that our 
approach to modeling, involving nonlinear indicial 
responses and generalized superposition integrals, could 
be made compatible with one based on nonlinear func- 
tional expansions. Here, we show that the form of the 
expansions is consistent with solutions of the unsteady 
Navler-Stokes equations and promotes the introduction of 
transform methods. 

Second, we discuss the issue of Frechet differen- 
tiability of the aerodynamic response. We have already 
noted that aerodynamic bifurcation is but one means of 
invalidating Frechet differentiability, and we have 
proposed that a theory for enumerating these means 
offers a possible alternative to bifurcation theory. 
Inasmuch as it would incorporate the latter, it would be 
a potentially more inclusive way of classifying steady 
and unsteady aerodynamic phenomena that are Important in 
flight dynamics applications. Here, we begin such a 
theory with an account of the means of invalidating 
Frechet differentiability that we have currently 
identified. 

Finally, we take up a potentially Important issue 
that emerges when we try to implement the accommodation 
within the mathematical model of bifurcation phenomena 
having time-dependent equilibrium states. Recognizing, 
for example, that with a Hopf bifurcation it is neces- 
sary to specify phase as well as amplitude of the 
periodic equilibrium state, one can frame the general 
question as follows: given that the past motion deter- 

mines the initial condition for the specification of an 
indicial response, how much of the past motion must be 
acknowledged in order to ensure that the equilibrium 
state of the response is adequately specified? As we 
shall see, time-dependent equilibrium states are never 
entirely free of dependence on the past motion, and the 
dependence becomes increasingly binding with each suc- 
cessive bifurcation. 


gravity v 0 (g) being an arbitrary function of time s. 
The angle of attack a is defined to be the angle 
between the resultant velocity vector and the chord line 
of the wing, so that a * tan^ffVgUJJ/Ug). It was 
argued that the response in lift at time t to a single 
pulse in a at time 5, ( 5 ^ < t) would have the form 

AC L (t) * AC L (t - 5 1 ,o(e 1 )) dir 

s 2 »„(t - €,)Ca(C,)J n AC, O) 

n 

whereas the response in lift at t to a pair of pulses 
in a at 5 , and 5 2 ( with 5, < 5 2 < fc ) would have the 
form 

AC L (t) = AC L (t - ^.eU,))^ * AC L (t - 5 2 ,«(5 2 >) dir 


* 4 V^int.Z 


(2) 


with 


A 2 C, 


'int,2 


^ **mn (t " e 2’*2 " C 1 > C 2> 1 
m,n 

x [0(5,)]" A5, A5 2 


(3) 


and it was observed that the process of adding pulses 
could be continued indefinitely in the same way. Going 
to the limit of a continuous distribution of pulses 
starting at time 5:0 led to a summation of multiple 
integrals having the form 


C. (t) * C, + C. + C. 

L L dir L lnt,2 L int,3 


c w • £ J *» (t ■ «. 

n o 

-V f [o(5 2 )] B d5 2 f 

‘“• 2 A 


(4) 


(5) 


b nn/ t ' 5 2 ,5 2 " 


x [a<« 1 )] n dC 1 


(6) 


NONLINEAR ALGEBRAIC FUNCTIONAL EXPANSIONS 

Reference (J) contains an analysis showing how a 
form for the nonlinear aerodynamic response in lift 
coefficient C^ to an arbitrary variation in angle of 
attack a could be constructed by means of a suitable 
summation of lift responses to pulses in a. To set the 
stage for what follows, the reasoning behind the form of 
the result obtained in (J) is summarized briefly below. 

Figure 1 illustrates the motion under study when it 
is referred to an X,Y coordinate system that is fixed 
in space. The wing moves away from the origin with 
fixed axial velocity U 0 (a slight change from the 
motion considered in (%) where resultant velocity was 
kept fixed), and at the same time is allowed to trans- 
late vertically with vertical velocity at the center of 


m,n,p o A 

i: 


[aU 1 )jPc mnp (t ’ 5 3' 5 3 ’ e 2'<2 ' « 1 )dt 1 

(7) 


Written as a nonlinear functional expansion, 

Eq. (4) represents the lift coefficient at time t in 
response to an arbitrary variation of angle of attack 
a over the time interval 0 to t. 
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The form of Eq. (4) could bo put on a firmer basis 
If It could bo shown to bo consistent with a solution of 
the Navier-Stokes equations. To that end, we shall 
write an analogous fora for the stream function *(x,t) 
governing the unsteady flow past a two-dimensional wing 
and inquire whether in fact it is consistent with a 
solution of the time-dependent Navier-Stokes equa- 
tions. Referring again to Fig. 1, let us transfer coor- 
dinates to the moving inertial frame (x,y) in which the 
fluid velocity components are respectively u and v, and 
where u aproaches the constant velocity U 0 as 
| x | * «. We confine the study to two-dimensional flow, 
which enables the Introduction of a stream function 
*(x,t). let us further require that the flow be incom- 
pressible and that the stream function satisfy the cor- 
responding Navier-Stokes equation, written in the con- 
venient form 


with 


= J - 5 1 ,*)o(5 1 )d5 1 


ft 

U| s I a 1 (t - (..xMc^dg. 

° y 

ft 

v, * -J a 1 (t - C 1 ,x)o(? 1 )dg 1 


(13) 


We see that with <i>j,Uj,Vf all having the form of con- 
volution integrals, we are able to introduce Laplace 
transforms to seek a solution of Eq. (12). Defining the 
Laplace transform as 


in ♦ uw + via - vv u = 0 
c x y 


with 


( 8 ) 


STfc[f(t)] = J e* st f(t)dt 


(14) 


“ = U y " V x * * ♦ 
U 1 


T y 

V = -e x 

Guided by the form of Eq. (4), we assign *(x,t) the 
form 


(where the superscript s and subscript t will be 
useful later to Identify the transform variable and 
(9) transformed variable in question) and letting 


♦( x , t) = e Q ( 


j) + Z J a n (t - «1* J 




n o 


Z f f 

m rt 


(•<«,)!" 


m,n o 


tf*[a(t) n ] = -*(s;n) 

^“[a^t.x)] = A^s;?) 
we get in Eq. (12) 

jW(s;1){sv 2 A, + u v 2 A. + v v 2 A, ♦ u A, 

1 0 1 O 1 O 1 

X y x y 


- n> 0 • vV Aj) s 0 


(15) 


( 16 ) 


y b (t - - e.,x)dc + ... (10) 

mn & * ' 

and ask whether the form is consistent with a possible 
solution of Eq. (8). After finding the forms of u, v, 
and u from Eq. (9) and substituting in Eq. (8), we 
attempt to match terms of like powers of a. Let us 
assume that prior to time zero the flow is steady and 
the angle of attack is zero. Then to O(a’) we have the 
initial flow + (x) which must satisfy the time- 
independent Nav?er-Stokes equation 

u o“o x + Vo y ■ v ’ 2 "o = 0 (11) 

To Ota 1 ) we have a linearized form of the time- 
dependent Navier-Stokes equation 

3 2 

rr id. * u w, + u,w + v w, + v.u - vV u, - 0 (12) 

«t l o 1 1 o o 1 I o 1 

“ * y / 


where, as required by Eq. (13), we have imposed the 
initial condition u^tO.x) * 0. For arbitrary a(t), 
the bracketed term must be zero, which provides a form 
for the solution of a^tt.x). 

To 0(o 2 ) we have 


it “2 + u o*2 + U 1“1 - u 2"o * v o“2 + V 1“1 

x x x y y 


V«w • W «- s 0 
£ O.. £ 


(17) 


with 


^ 2 r 4 

*2 3 J a 2 (t ‘ ♦ J o(5 2 )d5 2 

'o *o 

*f 2 <*(5 1 )b 11 (t - ? 2 ,i 2 - e 1 ,x)dt 1 


(18) 
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whereupon Eq. (23) becomes 


“2 = ’*2 


u 2 s *2 


V 2 = -♦2. 


In addition to terms similar in fora to those in 
Eq. ( 12), Eq. (17) contains terns consisting of double 
integrals and products of single integrals. Provided 
these can be shown to be transformable, Eq. (17) also 
will yield to Laplace transformation. The double- 
integral terms arise from the second contribution to * 2 
(Eq. (18)). Its Laplace transform takes the form 


r t t *22 1 = / 


f e -8t dt f o(e-)dt 2 f 
o *0 *0 




b 11 (t ‘ 5 2 ,5 2 ‘ *1 ,x)dC 1 


After reversals of the order of integration and the 
substitution t - c 2 3 u > 

*t [ *22 ] 3 f a(5 2 )e 352 d <2 f 2 •(«,)«, /” «‘ 3U 


* b 11 ( u te 2 - e 1 t x)du 


Letting 


*u [b 11 (u,t 2 * 5 1> X)I * B n (s,E 2 - £,;*) (22) 


we have 


**[♦ 22 ) 3 / a <«2 )e 2 de 2 °(«i) 

x B n (s,C 2 - ? 1 ;x)dt 1 (23) 

Recognizing that the inner integral is in the form of a 
convolution, we let 

3 1 t m ' s i w 

[B 11 (s,w;*)) = J e B^s.wjxJdw = B 11 (s,s 1 ;x) 


1 

^ t^22^ ~ 2*1 J — s 1 ;1Ms | ;1)B. 1 (i,s | ;i)ds 1 

»c-i« ‘ 

(26) 

Products of Integrals arise from terms such as, e.g., 

u.u, in Eq. (17), the Laplace transform of which has 
* ’x 

the form 

*t [u iWi 1 * f“ e" 3t dt^j" 6 a, (t - e 1 ,x)o(C 1 )d5 1 ) 

* [j* (t - c 2 ,5)a(? 2 )dC 2 J (27) 

Terms such as this can be handled by the rule (4) 

f,c t )p 2 <. - < 28 ) 

Using Eq. (28) and noting that each of the product terms 
is itself a convolution, we get for the transform of 

“’“’a 

1 /*c*i- 

^(u,-, 1 1 2vl J t ^(Sj, 1 Ms - s 1 , 1 )A 1 (s,;x) 


V A 1 (s - s 1 ;x)ds 1 


So Eq. (17) yields to Laplace transformation. We get 

^(s; 2){3V 2 A. * u V 2 A * v V 2 A + » A_ - w A 
c x y x y y x 


- vV Aj} ♦ 


_1_ f c+1 
2*1 f 


^(s^lWls - SjjIMsv B 11 (s,a 1 ;x) 


* vV ♦ vV * “o B n - “o B 11, ■ vv4b ’ 1 

a y a y y * 

♦ A. (s 1 )V 2 A 1 (s - s ^ ) - A^ (s 1 )V 2 A 1 (s - s^Jds, s 0 

y x x y 


so that 


rh 

J o(c 1 )B 11 (s, 5 2 - 5 1 ;*)d« 1 


. « S.£_ 

2 ii J • rf (*^i1)B 11 (s,s 1 ;x)e ds^ (25) 


Again, with o(t) arbitrary, each of the bracketed 
terms mist be Independently zero, which provides forms 
for the solution of a 2 (t,x) and b^tt^ ,t 2> x). We see 
that a 2 (t,x), and in fact, each subsequent a i (t,x), is 
governed by the same linear fora as that of a 1 (t,x) . 
Each will have the same set of elgensolutlons, but their 
combinations will differ. The homogeneous operator on 
b j y also has the same fora (as will that on each subse- 
quent b ij\ but bjj depends additionally on a^(t,x). 

It is clear that the procedure can be continued to 
higher powers of a. At each stage, the multiple 
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integrals can be transformed to functions of a single 
transform variable and forms will be available for the 
solution of the various transformed coefficients. Thus, 
although the problem of satisfying boundary conditions 
remains to be addressed. It appears that the fora of 
Eq. (10) will be consistent with a possible solution of 
the Navier-Stokes equations. It Is in effect a power- 
series expansion In a around a = 0, and should be 
valid for at least a small range of a In the vicinity 
of a s 0. The ways in which It can fail brings up the 
subject of Frechet differentiability. 

FRfiCHET DIFFERENTIABILITY 


To study the Issue of Frechet differentiability, 
let us consider the general problem of forming the indi- 
clal response Imagining that it can be obtained from 
solutions of the Navier-Stokes equations. We need to 
consider two motions (cf. (J) and Fig. 2). In the 
first, the body undergoes the motion under study a(c) 
from time zero to time 5 - t. Subsequent to x,a is 
held constant at o ( t ). The motion history is pre- 
scribed as <^(5): 


c^U) = a(0 ; 0 < C < x 

s «(t) ; 5 2 t 


(31) 


We ask for the response in, say, surface pressure 
p(x 3 ,t), which can be written in functional notation as 

p(x s ,t) s pt«» 1 (C) ;* 3 J = p[«(t);t»T,x 3 ] (32) 

Let us assume that a solution for p(x 3 ,t) can be 
found from the appropriate two- or three-dimensional 
time-dependent Navler-Stokes equations. In the second 
motion, the body undergoes the same angle-of -attack 
history a(c) up to time t . Subsequent to t , the 
angle of attack is again held constant, but is given an 
incremental step change ia over its previous value 
of o(x). Following (2), we represent the second motion 
as 


with 


a 2 U) — ajU) + en 


(33) 


e = ha 

n = 0 ; 0 < t < x 

= 1 ; tS( 


(34) 


The surface-pressure response to this motion is 

p(x 3 ,t) = pfajU) + *n;x s 3 (35) 

Again we assume that a solution can be found. The indl- 
cial response in surface pressure is the Frechet deriva- 
tive, which is defined as 


Ap(*,.t) 

11m r 2 = lim 


pt^U) + cn;x 3 J - p[a 1 (c);z g ]i 


Aa+0 


La 


€♦0 


= 57 pc°iU> * e ";* s i £I o 


= p'(a 1 (e);x 3 ]n i P a t<»U);t,T,x 3 ] (36) 

Let us assume that the derivative exists for every a 
that occurs in a motion over some time Interval. We 
obtain the surface pressure response to the angle-of- 
attack variation by forming a generalized superposition 
Integral : 


p(x 3 ,t) — p(t,o(0);x 3 ) + j p a (a(E);t,x,x 3 ] ^ dx (37) 

Validity of the approach thus rests on the Frechet dif- 
ferentiability of the surface-pressure response over the 
Interval 0 < x < t. 

As we have noted (2), a theory for enumerating the 
means of invalidating Frechet differentiability of the 
aerodynamic response offers the possibility of a more 
inclusive way than bifurcation theory is of classifying 
steady and unsteady aerodynamic phenomena that are 
important in flight-dynamics applications. At present 
we have identified three distinct means. 

Aerodynamic Bifurcation 

We define aerodynamic bifurcation as the replace- 
ment of an unstable equilibrium flow by a new stable 
equilibrium flow at a critical value of a parameter. 

The onset of instability of the equilibrium state is a 
linear problem, capturable, at least in principle, 
through analysis of a linearized form of the time- 
dependent Navier-Stokes equations. For two-dimensional 
flow, the governing equation for the perturbation stream 
function is the same as that of Eq. (12), with the 
exception that the basic state (u 0 ,v Q ,u 0 in Eq. (12)) 
must be the equilibrium state that is existent at the 
angle of attack a(x) in question. For this reason, 
analysis based on the solution form Eq. (10) will not 
reach far enough to identify bifurcation phenomena with- 
out further generalization. 

An example that illustrates aerodynamic bifurcation 
in a context relevant to current interests is the slen- 
der body of revolution at high angle of attack. We 
shall describe this problem heuristically by adapting 
the impulsive flow analogy, both to demonstrate the 
continued usefulness of the analogy as a conceptual 
device and to make a point that will be important 
later. The analogy has been used extensively to help 
explain the presence of the steady synmetric vortices 
that are found on the leeward of slender wings and 
bodies at high angles of attack. To review briefly, we 
refer to Fig. 3 where now the moving axis system (x,r,8) 
is fixed to the body nose and the body is sinking with 
uniform velocity v 0 . Fixing attention to a vertical 
plane fixed in space, e.g., the plane X = 0, and ignor- 
ing the growth of the body, one observes the flow about 
a cylinder that appears to have started from rest and is 
sinking at a uniform rate v c . On the assumption that 
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the flow In the vertical plane la essentially two- 
dimensional, one is able to relate the steady growth of 
the vortices In the body-axis system with distance x 
along the body to their time-dependent growth behind the 
impulsively started cylinder, as viewed in the X s 0 
plane. The latter problem has been studied extensively 
both experimentally (£) and via numerical computation 
(6) . The first appearance of separation on the leeward 
ray, followed by the growth of the primary vortices, 
then the secondary vortices, and so forth, has been 
thoroughly documented. We also know that beyond a crit- 
ical Reynolds number (Re ■ 50 based on v c and on body 
diameter), the equilibrium flow becomes time-dependent, 
leading to the periodic shedding of vortices in the 
wake. This phenomenon also has been studied exten- 
sively. The work of Nlshloka and Sato (7) is particu- 
larly noteworthy, inasmuch as it demonstrates clearly 
that the onset of periodicity occurs at a critical 
Reynolds number and constitutes, in effect, a Hopf 
bifurcation. Ue expect, then, that at a critical 
Reynolds number a perturbation velocity within a fixed 
plane will have behavior that can be expressed Initially 
as 

v * g(r,9)e l “ fc (38) 

where if v is radial velocity, then g(r,8) is anti- 
symmetric with respect to e = 0. Now, we have assumed 
that the axial flow induced by the cylinder as it moves 
through the plane is more or less Independent of X. As 
we shall note again later, a nominally two-dimensional 
flow will permit a wave-form perturbation solution in 
the direction of the third coordinate. A more general 
form of the perturbation velocity v is then 

v = g(r,8)e i “ t e UX (39) 

Transforming this solution to the body axis system 
with X = x - U Q t, we get 

iX[x-U t+(»/X)t] 

v = g(r,e)e (40) 

and we see that a stationary perturbation solution is 
possible in the body-axis system when the wave speed 
u/X equals the axial speed U Q . Experimental results 
indicate that a stationary asymmetric flow is the pre- 
ferred solution over a considerable range of angle of 
attack. He conjecture that it is easy for an equality 
between U Q and the wave speed ui/x to be achieved for 
the following reason. Owing to the presence of the 
vortices, the axial velocity profile has a considerable 
overshoot beyond the free-stream value U Q . The wave 
speed u/X must lie between the extremes of the axial 
velocities available and hence can easily match U Q . On 
the other hand, as o ♦ 90°, X ♦ 0, a match becomes 
Impossible, and the perturbation form reverts to 
Eq. (38) in either axis system. 

The analysis leads to the following conclusions. 

(1) Having demonstrated its ability to fit observations, 
the impulsive flow analogy remains a useful conceptual 
device for understanding complex three-dimensional 
flows. (2) Accepting (1) encourages the belief that 
Instabilities in complex three-dimensional flows may be 
essentially local phenomena which will yield to 


localized analyses. (3) The form of bifurcation, and 
hence its classification, may depend on the coordinate 
system in which it is observed. Frechet differentia- 
bility on the other hand, remains invariant across 
coordinate transformations. 

Loss of Analytic Dependence on a Parameter 

The following refers principally to time-invariant 
equilibrium flows. The variation with a parameter 
(e.g., roll angle *) of the aerodynamic response may 
develop a fold at a critical value of the parameter * 0 , 
so that the slope of the response becomes infinite 
there. A Jump in the response necessarily ensues to 
another branch of the folded curve with an infinitesimal 
increase in 8 beyond 8 0 , and hysteresis follows on 
the return route. 

Folds need not be the result of bifurcation; no new 
branches of equilibrium flow solutions need appear. On 
the other hand, folds may be the indirect result of 
bifurcation. An example based on the preceding study is 
the slender delta wing. Just as for the slender body of 
revolution, bifurcation to an asymmetric flow pattern 
occurs at a critical value of a, beyond which a finite 
value of the rolling-moment coefficient C t may exist 
at zero roll angle. As shown in Fig. 4, both a positive 
and a negative value of C t is possible, depending on 
the sense of the asynmetry. The solution C t = 0, cor- 
responding to the unstable symmetric flow, also exists 
at 8-0, but is itself unstable. So folds must occur 
to accommodate the three solution points. 

Disconnected Bifurcation 

By disconnected bifurcation we mean the existence 
of an isolated branch of equilibrium solutions that is 
not connected to other solution branches at bifurcation 
points. The sources of this form of invalidation are 
still obscure although evidence abounds of its existence 
in experimental studies. It is clear that the state of 
the flow at the origin of the indlcial response e s r 
must be involved. Disturbances there must be large 
enough and of an appropriate form so as to divert the 
response toward a new attractor, that Is, toward an 
alternative branch of equilibrium solutions. 

As we have already seen, a source to which nomi- 
nally two-dimensional flows should be particularly sus- 
ceptible is the presence of three-dimensional distur- 
bances in cellular patterns In the direction of the 
third coordinate. As candidates for study in this light 
we cite the following two possibilities. The first one 
concerns flow at the stagnation line of a cylindrical 
body. The presence of a regular pattern of disturbances 
along, e.g., the leading edge of a straight wing has 
been known for a long time. Attempts to explain the 
pattern's presence on the basis of stability theory have 
been neither consistent nor convincing. In a comprehen- 
sive and useful review of the subject, Morkovin (8) 
concludes that there is still no viable explanation of 
the phenomenon. The second possibility could be axisym- 
metric vortex breakdown. 3 A solution here on the basis 


3 Joint work with D. Weihs, Technion, Israel. 
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of a disconnected bifurcation would represent an exten- 
sion to viscous flows of Brooke Benjamin's notion (J) of 
"conjugate flows" that he Introduced as an aid to the 
explanation of inviscid axisynrntrlc vortex breakdown. 

We believe that an extension of our use of nonlinear 
functional expansions that accounts for cellular pertur- 
bations In the third coordinate direction will provide a 
focus on the sources that may be available to sustain 
such perturbations in the Initial flow. 

MEMORY EFFECTS 

We have postulated six major subdivisions in the 
form of the aerodynamic response by means of a set of 
sketches (J) , and these are reproduced In Fig. 5. Of 
particular interest here are Figs. 5(d) and (e), deplet- 
ing time-dependent equilibrium states. Hopf bifurcation 
Is indicated in Fig. 5(d), wherein a formerly stable 
time- Invariant equilibrium state Is replaced by a time- 
varying periodic equilibrium state. The equilibrium 
state resulting from a Hopf bifurcation is often suc- 
ceeded by bifurcation to a quasi-periodic equilibrium 
state with further Increase of the relevant parameter 
(here, angle of attack). The quasi-periodic state may 
be succeeded by an aperiodic (chaotic) equilibrium 
state. 

We Indicated (£) how these bifurcations could be 
accommodated within our mathematical model. In the 
implementation, however, an important additional issue 
first arises with a Hopf bifurcation when we recognize 
that it is necessary to specify phase as well as ampli- 
tude in order to completely determine the periodic equi- 
librium state. The general problem can be posed as 
follows: given that the past motion determines the 

initial condition for the specification of an lndiclal 
response, how much of the past motion must be acknowl- 
edged in order to ensure that the equilibrium state of 
the response is adequately specified? We see that there 
will be a definite sequence, in which at least one piece 
of additional information about the past motion must be 
supplied with each successive bifurcation. (The 
sequence is illustrated in Fig. 6.) Let us note, Inci- 
dentally, that with a- linear system, no information 
about the past motion is required in the specification 
of an indlelal response. This is, in effect, the most 
general definition of a linear system. For a nonlinear 
system in which the equilibrium state is time-invariant, 
the level of the past motion at the origin of the step 
is all that is required to specify the equilibrium state 
(Fig. 6(a)). For the time-periodic equilibrium state 
that replaces the time-invariant state with a Hopf 
bifurcation, the amplitude and frequency of the periodic 
state are again determined by the level of the past 
motion, but an additional piece of information, the rate 


of the past motion at the origin of the step must be 
supplied to fix the phase of the equilibrium state 
(Fig. 6(b)). If the next bifurcation results in a 
quasi-periodic equilibrium state, another derivative of 
the motion at the origin of the step will be required to 
specify the additional phase (Fig. 6(c)). It is clear 
that a chaotic equilibrium state will require informa- 
tion about the complete past motion to ensure a complete 
specification. This is how "sensitive dependence on the 
initial conditions," the signature of a chaotic state, 
will manifest itself in the modeling of the tndlcial 
response. We recognize the novel appearance of a turbu- 
lence modeling question: Is it possible to model the 

past motion so as to suppress the chaotic part of an 
equilibrium state? The necessity of coping with this 
problem in the context of flight dynamics studies demon- 
strates once again the convergence of interests that is 
occurring among workers in heretofore widely separate 
fields. 
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Fig. 1 Wing maneuver referred to space-fixed (X,f) and 
moving (x,y) coordinates. 
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Fig. 2 Formation of lndiclal response in surface 
pressure. 
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(c) QUASIPERIODIC EQUILIBRIUM 
STATE 


Id) CHAOTIC EQUILIBRIUM 
STATE 


Fig. 6 Memory effects. 
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